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NOMENCLATURE * 

A Source strength coefficient 

b Slope of velocity profile in shear layer, du/dz or dU/dr 

B (b/*)K;[ in Cartesian coordinates, or 

(b/w)K in cylindrical coordinates 
c Speed of sound 

f First solution of ordinary differential equation in 

shear layer 

g Second (independent) solution 

i V-l 

k(z) Adiabatic constant for given z 

wave number in x-direction (Cartesian coordinates) 

k^ Wave number in y-direction (Cartesian coordinates) 

k Wave number = Vk^^+k^^ (Cartesian coordinates); 

Wave number in axial direction (cylindrical coordinates) 

K Reduced wave number in x-direction = k /« f (Cartesian 

1 coordinates) 1 

K 2 Reduced wave number in y-direction = (Cartesian 

coordinates) 

K Reduced wave number = k/w’ = Vk ^+K ^ (Cartesian coordinates) 

12 

Reduced wave number in axial direction = k/w 1 (cylindrical 
coordinates) 

M Mach number (based on jet temperature unless subscript a 

appears) 

m an integer 

n an integer 

p Pressure 

*The secondary symbols not appearing in this list are defined 
where they are introduced in the analysis. 
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NOMENCLATURE * 

A Source strength coefficient 

b Slope of velocity profile in shear layer, du/dz or du/dr 

B (b/“')K 1 in Cartesian coordinates, or 

(b/o>)K in cylindrical coordinates 
c Speed of sound 

f First solution of ordinary differential equation in 

shear layer 

g Second (independent) solution 

* s 

i V-l 

k(z) Adiabatic constant for given z 

kj Nave number in x-direction (Cartesian coordinates) 

k 2 Wave number in y-direction (Cartesian coordinates) 

k Wave number = (kj^+kg 2 (Cartesian coordinates);. 

Wave number in axial direction (cylindrical coordinates) 

K Reduced wave number in x-direction = k /u>' (Cartesian 
1 coordinates) 1 

K 2 Reduced wave number in y-direction = k 2 / <*> * (Cartesian 
° coordinates) 

K . Reduced wave number = k/w’ = i K 2 +K 2 (Cartesian coordinates) 

JL Z 

Reduced wave number in axial direction = k / u' (cylindrical 
coordinates) 

M Mach number (based on jet temperature unless subscript a 

appears) 

m an integer 

n an integer 

p Pressure 

*The secondary symbols not appearing in this list are defined 
where they are introduced in the analysis. 
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AP 


Pressure increment from acoustic waves 


Ap^ Mean-square pressure from acoustic waves 
r Radial cylindrical coordinate 

R Radial spherical coordinate = x 2 ^ + ^ +z ^ - >. 

R* Hyperbolic radius = / x 2 +( l-M s 2 a T ) r 2 

R Gas constant 

R.P. Denotes "real part of" 

s Condensation 

Sj Strouhal number 

T Absolute temperature 

t Time 

u Axial velocity 

l 

U(z) or U(r) defines jet velocity profile 
v 'Lateral or tangential velocity 
w Normal or radial velocity 

x Streamwise coordinate (coordinates fixed in source) 

X 2 Streamwise coordinate (coordinates fixed in ambient air) 

y Lateral coordinate 

z Normal coordinate 

a^, Ratio of jet temperature to ambient temperature 

y Ratio of specific heats 

8 Half-width of uniform velocity region 

$ co'(z-z s ) when z - . z Q ; ^'(z-Zq) when z g r. z Q 

v (r- r s ) when r g :> r Q ; uP (r-r Q ) when r g . r Q 

6 Angular position in far field measured from upstream jet 
axis 

p Fluid density 
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0 Velocity potential 


Angle about x-axis (cylindrical coordinates) 
Generating frequency 

w/c 

Subscripts 

a refers to ambient 

c refers to critical 

j refers to jet 

s refers to source 
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1. INTRODUCTION 


Objectives 

This analysis is not intended as a comprehensive theory 
for the prediction of noise produced by turbulence in jets. 
Instead we concentrate on one aspect of this complex problem, 
the transmission of acoustical disturbances from the interior 
of the jet through a shear layer and into the far field.* 

These acoustical disturbances are generated by mathematically 
defined point sources drifting with the local fluid. We ne- 
glect temporarily the more difficult problem of identifying 
these mathematical sources with the turbulence which they are 
intended to represent. 

The problem we have formulated can be treated with some 
precision, and with comparatively few assumptions. In spite 
of these simplifications the results obtained are often not 
predictable by simple intuitive ideas. This is in fact the 
main advantage of such an approach, that, for those attempting 
to construct more comprehensive theories of jet noise, it 
should provide facts not intuitively evident. 

Description of Model 

The noise generator chosen is a sequence of transient 
sources drifting with the local fluid. The jet (either two- 
dimensional or circular cylindrical) extends to infinity up- 

*This includes both "convection" and "refraction" effects 
which are in general inseparable. 
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stream and downstream, with velocity profile independent of 
streamwise position. Thus the large velocity gradients across 
the jet are accounted for, and the smaller gradients in the 
streamwise direction are neglected. It seems reasonable to 
expect that the major refraction effects will be shown by such 
a model. 

In making this analysis we consider that turbulence ("self- 
noise") is the only true originator of noise, and that "shear 
noise", being composed of linear terms, is a part of the trans- 
mission process. 

It must be emphasized that this is not a stability anal- 
ysis. We deal with a distribution of turbulence in the jet 
which is essentially independent of time. This steady-state 
situation is the end result of the action of instabilities. 

We then attempt to find the effect of the jet mean velo- 
city profile on the transmission of acoustic disturbances from 
one element of turbulence through the jet and into the far 
field. The scattering effect of other elements of turbulence 
is neglected. 

Although supersonic as well as subsonic jets can be 
treated by the methods discussed in this report, the supersonic 
cases present some anomalous features which require more careful 
study. The results given here are therefore confined to sub- 
sonic jets. 


» 
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2. DESCRIPTION OF METHODS 
Generality of Methods * 

All of the tasks outlined in this report have in common 
certain necessary steps and concepts. The noise may be gener- 
ated by sources (i.e. monopole sources), dipoles or quadrupoles. 
The noise generators may be on center, off center in a uniform 
velocity region or, in some cases, in the shear layer. Temper- 
atures may be ambient or varying across the jet. The jet itself 
may be two-dimensional or circular. Still, with only minor mod- 
ifications certain basic ideas apply, and we review them briefly. 

For convenience we will describe the jet as being two- 
dimensional. The velocity profile must be independent of 
streamwise and lateral positions. For simplicity we will say 
"source", and speak of a single (constant velocity gradient) 
shear layer on each side of it, but other singularities can be 
used, and more complicated shear layers can be treated. 

The Source in a Jet 

It is first assumed that the pulsating source is at rest 
in a completely stationary homogeneous fluid of infinite ex- 
tent. The origin of coordinates is fixed in the source, and 
the conventional wave equation applies. The velocity potential 
produced by the source is readily expressed as the double inte- 
gral of the-velocity potentials for all reduced wave numbers, 
and in the x (streamwise) and y (lateral) directions. 

*See references 1-5- 
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This corresponds to the decomposition of the source potential 
into an infinite set of plane waves (and exponential disturb- 
ances) . 

The source and some portion of the surrounding fluid at 
rest must next be confined within a jet. From coordinates 
fixed in the source this jet is seen as shear layers flowing 
past on each side, and outside the shear layers the ambient 
fluid flowing by at a fixed velocity. (See Fig. 1.) To ac- 
complish this insertion of the source into the jet one must add 
to its original velocity potential the potentials for upward- 
moving and downward-moving waves reflected off the shear layers. 
These two reflected wave amplitudes are as yet unknown. 

The Shear Layers and Ambient Air 

In the shear layers the conventional wave equation does 
not apply. The correct partial differential equation is de- 
rived, and by assuming periodic solutions in the streamwise and 
lateral directions an ordinary differential equation is obtained 
in the coordinate normal to the shear layers. The ordinary 
differential equation is solved by power series expansions 
about the singular point, and two independent solutions are ob- 
tained (in each layer) with amplitudes as yet unknown. 

In the ambient air the conventional wave equation applies 
again (for coordinates fixed in the ambient fluid). Only out- 
ward moving waves need be considered, so (in each ambient re- 
gion) one solution of unknown amplitude appears. 

There are now eight unknown amplitudes to be determined 
and four boundaries between fluid layers. Across each of these 
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bounaaries pressure and displacement must be continuous, yield- 
ing the necessary eight equations to make the system determinate. 
A double integral (over and K^) for velocity potential in the 
ambient air is then obtained as a function of source strength 
and frequency. 

The Far Field 

Only in the far field can this double integral be evalu- 
ated easily. There the integrand generally consists of a slow- 
ly varying function multiplying a function which oscillates 
rapidly about zero. Such functions interact weakly, and sig- 
nificant contributions to the integral occur only when (a) the 
slowly varying function becomes rapidly varying (i.e. in the 
neighborhood of singular points) or (b) when the rapidly oscil- 
lating function ceases to be rapidly oscillating (i.e. in the 
neighborhood of stationary phase points). Evaluation of the 
integral yields the pressure in the far field for coordinates 
fixed in the source. 

Transient Sources and Retarded Coordinates 

A more useful result would be the far-field mean-square 
pressure produced by sources in a localized region, say immed- 
iately behind a jet nozzle. To obtain this we consider a se- 
quence of transient sources, each originating at the same point 
relative to the ambient air or nozzle. As one source disappears 
after drifting downstream with the fluid, a new one (with random 
phase relative to the first) appears at the upstream point. The 
mathematical analysis of this process (kef. 3) is rather tedious. 
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However the practical application conforms to a simple rule if 
it is assumed (as in the present analysis) that the transient 
sources have a lifetime of many cycles. The rule is that the 
mean-square pressure in source coordinates should first be 
formally transformed to retarded coordinates, fixed in the 
nozzle. (See Fig. 2.) The result must then be multiplied by 
• 1 + MX 2 /R 1 t where M is the source Mach number relative to the 
ambient air, X 2 is the streamwise distance and R the radius to 
the far-field observation point.* (The rule holds also for the 
supersonic case if the |1 + Mx^/Rj factors are enclosed by ab- 
solute magnitude signs.) 


Contrast with Lighthill's Procedure 

It may be useful to point out a few of the features in 
which our procedure differs from that of Lighthill^. 

In our analysis only the non-linear terms (the turbulence 
itself) are assumed to generate noise. All linear terms, in- 
cluding those commonly associated with "shear noise", are in- 
cluded in the transmission process. 

Also, our sources are at rest in the immediately sur- 
rounding fluid and so are defined unambiguously. Lighthill's 


* The factor 


1 + MX 2 /R is often written 1 + M cos0 where 
9 is the angular position of the far-field observation point 
relative to the jet axis, 9 = 0 being measured upstream. 


0 = 7 t downstream . 

The travel distance of each transient source and the thickness 
of the jet are of course negligible compared to R. 
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sources are in motion relative to the neighboring fluid, and 
belong to one of the (in principle) infinitely many classes 
of moving sources 5 . 

In Lighthill's model of the noise-producing jet no fluid 
moves. Instead, an array of acoustical point sources passes 
through the fluid, all sources moving at precisely the same 
velocity. It is this simplification of the model which leads 
to the elimination of source and dipole noise, and leaves 
quadrupoles as the dominant noise generators in Lighthill's 
analysis. 
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3. DISCUSSION OF DIPOLES AND QUADRUPOLES 

Since quadrupoles are frequently assumed to be the pri- 
mary noise generators in a jet, our analyses have included 
such singularities located on the center-line of a two- 
dimensional jet (see Appendix). 

However, in examining the construction of dipoles and 
quadrupoles in a shear layer we were led to the conclusion 
that there is no obvious reason for excluding the simple and 
highly efficient monopole source as a primary noise generator. 

A discussion of these ideas follows. 

We visualize a dipole as being the limit approached when 
a source and sink* of equal strength are brought together, the 
strength being inversely proportional to the separation dis- 
tance. In a region of uniform fluid velocity (such as assumed 
at the centerline) this presents no conceptual problem. The 
source and sink, each drifting with the fluid, travel at the 
same velocity and maintain the same relative positions. In a 
shear layer this may not be true. Roth the length and direc- 
tion of the line connecting source and sink may vary. One 
might suppose that in the limit, when source and sink coincide, 
this would be unimportant. However for a constant velocity 
gradient in the z-direction the percentage change in length 
and the rotation of the axis are independent of the z displace- 
ment, and may be 40 percent and 90 degrees, respectively, when 
source and sink travel one jet diameter (or less). Another, 
and possibly more important, problem is that the faster moving 

* 


a pulsating source 180-degrees out of phase 
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source generally produces much greater pressures in the far 
field. As the source and sink approach each other this effect 
might be lost except that their strengths are being increased 
in the process of dipole construction. To illustrate this 
possibility we consider d(0 s )/dz as compared to d(z0 s )/dz, 
where 0 S is source potential and z0 g crudely weights the effec- 
tive strength for radiation. The former is proportional to the 

potential of an ordinary dipole, but the latter contains a 
residual source potential. 

What we now suggest is that the absence of any net source 
strength in the entire jet does not necessarily result in 
purely dipole radiation, nor does the absence of a net dipole 
strength result in purely quadrupole radiation in our analysis. 
Thus one of the commonly used arguments in the exclusion of 
source and dipole radiation apparently does not apply here. 
There may be other and better reasons for choosing quadrupoles 
in some instances. Until these are clarified, however, we 
propose to deal with the simpler source noise generation in the 
subsequent tasks of this study. 
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4. THE SOURCE IN THE SHEAR LAYER OP 
A TWO-DIMENS IONAL JET 

A. DEVELOPMENT OF EQUATIONS 

Symmetric and Anti-Symmetric Solutions 

To analyze the case of a source located in the shear layer 
of a jet (Fig. 3), we first assume a narrow region (width 26) 
of uniform velocity (U g ) within the shear layer. The source is 
located in this uniform-velocity region, and is traveling with 
it, so that the source is at rest relative to the immediately 
surrounding fluid. In a coordinate system moving with the 
source the jet velocity distribution is then as sketched in Fig. 
4. This problem can be solved by the methods briefly outlined 
in Section 2. The thickness of the uniform-flow region is then 
allowed to approach zero. 

In order to reduce the number of simultaneous equations 
to be solved we utilize in practice the properties of symmetri- 
cal and anti-symmetrical solutions. First, we assume two sources 
symmetrically disposed about the centerline of the jet and pulsa- 
ting exactly in phase. This problem is completely symmetrical 
and can be solved in the half- jet by requiring no displacement 
at the centerline. Next, the anti-symmetrical problem is solved 
considering the symmetrically disposed sources to be pulsating 

180-degrees out of phase. This problem is also solved in the 

t 

half- jet, but using zero pressure at the jet centerline. The 


superposition of velocity potentials for these two cases gives 
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Fig. 3 TWO-DIMENSIONAL JET WITH 
SOURCE IN SHEAR LAYER 
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x- = x-U t 
£ s 




z = z 



Velocity and distance are 
positive to the right. 


HALF-JET 

L_ 

Region 2 


c 

Region IB 

) Source Region 6 

Region 1A 
Region 0 


Jet Centerline: Symmetric case - zero displacement 

Anti-symmetric case - zero pressure 


Fig. 4 JET GEOMETRY AND VELOCITY PROFILE 


(Coordinates Fixed in Source) 
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the result for a single source (of double strength). The mean- 
square pressure is then found for a single source of the proper 
strength. 

Since the jet is symmetrical we can expect, on the average, 
to find a similar source on the opposite side of the jet, but 
with random phase relative to the first source, and we add their 
mean-square pressures. (Note that the earlier use of in-phase 
and out-of-phase source pairs was a purely mathematical exped- 
ient for constructing a single source on one side of the jet.) 


Coordinate Systems 

The basic coordinate system (x,y,z) is fixed in the source, 
with x the streamwise, y the lateral and z the normal coordinate. 
Coordinates fixed in the ambient air are X2 ,y,z* and those fixed 
in the uniform flow region at the jet center are Xq,v,z, where 

x 0 = x-U_t 

2 S (1) 
x 0 = x+(Uj-u s )t 

Here U„ and U ; are, respectively, the absolute magnitudes of 
the source and jet velocities relative to the ambient air, and 
t is time. 


The Central Uniform-Flow Region 

In Region 0 (see Fig. 4) the flow is uniform, and the con- 
ventional wave equation is satisfied in coordinates fixed in the 
f luid: 


0 

* 0*0 






0 


( 2 ) 
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where 0 is the velocity potential, c is the speed of sound and 
t is time. 

To satisfy the wave equation, provide for upward- and 
downward-moving waves and satisfy the condition of zero dis- 
placement (or zero pressure) at the jet centerline we choose 0 
to be 


0 Q = Aw' R.P. R 


exp 


iwUK,x+Kj-ct+zm+K,(M.-M )| 

12 L 1 J s - 1 J 


2 v 2' 


expj^i a,' (^x+^y-ct-z ^[l+KjCMj-Mg)] 2 -K 2 ) 


(3) 


where x = x -(U.-U )t. As mentioned earlier we are dealing with 
0 j s 

two cases here. The symmetric case represents a source at z = z g 
and a source at z = -z , and in the half- jet at z = 0 there is 
required 

(d0 o /dz) = 0 (4) 

The anti-symmetric case represents a source at z = z g and a 
"sink" at z = -z g , and in the half jet at z = 0 


(d0 Q /d t) 


XQ=const . 


0 


(5)* 


This latter case is indicated by the minus sign in parentheses 
before the second term of Eq. (3). 

2 

In Eq. (3) , A = (source strength)/8 tt 9 the source strength 
being the maximum volume introduced per unit time. w T = w/c 
where u> is the generating frequency in radians per second and 

^pressure Cap) = - p 0 Ot = - P (d *V at) x 0 =const. 
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c is the speed of sound. R is a complex reflection coefficient, 
R = R'+iR". K x = kj/u)', K 2 = , K 2 = K * + K^ 2 . k and k 2 

are respectively the wave numbers in the x and y directions. Mj 

is the maximum jet Mach number (M . = U./c) and M is the source 

J / s 

Mach number (M = U /c). R.P. denotes "the real part of". (The 
s s 

jet is here assumed to be at ambient temperature, so that 
p and c retain constant values throughout the flow field.) 


The Source Region and the Ambient Air 
In Region 8 (see Fig. 4) we have 


0 = Aw* R.P.<! 


expiiw' (K x+K y-ct+ z-z 
! 1 2 s 


/] L-r 2' ) j /yl-K 2 ' 


i «>« (K x+K y-ct-z /l-K 2 ) | / ti-K 2 
x ^ , 


+ P exp 

t. X 6 _! 

+ J expjji »* (K 1 x+K 2 y-ct+zVl-K 2 )J /vl-K 2 L 


( 6 ) 


where the first term, integrated from - M to +oo in both K^ and 
K 2 t gives the potential for a source in a homogeneous fluid of 
infinite extent. The second and third terms account for down- 
ward and upward moving waves reflected off the boundaries of the 
8-region. P and J are respectively P‘+iP" and J’+iJ", complex 
reflection coefficients. 

In region 2 (Fig. 4), the ambient air, we have 


0 2 = Am' R.P. S exp 


a u>’ 


Kp (x 2 +b s t ) +K 2 y- c t i ( z-z x ) V( 1 -K x M s > 2 -K 2 J > 

(7) 


where the lower sign is used for reversed waves (c-K U < 0 and 

x s 

2 2 

(1-K 1 M ) > K ). Otherwise the, upper sign is used (ordinary 

X . s 

7 7 

waves when c-K^Ug > 0 and (1-K^Mg) > K ; exponential decay when 
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2 2 — 

(1-K^Mg) < K ). S is a complex transmission coefficient here 

defined as S = S'+iS" for convenience in satisfying boundary 
conditions in complex form. (In some of our previous analyses 
S was defined as S M -iS' , but this change causes no confusion in 
the final result where only S'^+S"^ appears.) 


The Shear Layer Regions 

So far we have considered only the regions of uniform flow. 
The shear layer regions must be treated differently, using s, 
the condensation, and w, the vertical velocity, instead of the 
velocity potential 0. The conventional wave equation is re- 
placed by a more complicated form^ which, upon the assumption of 
simple harmonic forms in the x and y directions, yields an or- 
dinary differential equation (Eq. (17), given later) in the z 
(or C) direction. 

In Region IB (see Fig. 4) 


s 

w 


R.P. Qg( Cg) exp £i w’ (KjX+K^y-ct )j 

R.P. [-ic/(l-Btg) Qg exp[i< u »(K 1 x+K 2 y-ct)] 


( 8 ) 

(9) 


where 

B = (b/oi)K 1 and b = dU/dz (10) 

U is the local shear velocity in the layer. £ is a coordinate 

B 

in the z direction defined by 


C B = w' [z - (z s +8) 


( 11 ) 


Q is a complex function of C and Q 
B B Dt. 


is the derivative of this 
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function with respect to C . (Note that R, P, J are complex 

B 

coefficients but independent of the z (or £) coordinate.) 

Qg is composed of solutions of the ordinary differential 
equation which applies in a shear layer (Eq. (17) ) as will be 
indicated later. 

In Region 1A, similarly. 


s = R.p. Q a U a ) exp[i w’ (KjX+Kgy-ct). (12) 

w = R.P. | -ic/(l-BC )• Q, exp[i w , (K 1 x+K y-ctj"! (13) 

^ A : A^ JL. c* 

where £ A is a coordinate in the z direction defined by 

C A = ««[z - ( 2 s - 8 )l (14) 

Now, if 8—0, C A — C , where 

C = ®'(z-z s ) (15) 

We may write 


Q A = *A f(C) + b A g(?) 

% = a B f(C) + b B g(0 


(16) 


where f(C) and g(C) are solutions of the shear layer differen 
tial equation* 


(1-B£)F CC + 2BF^ + (1-B£)| (1-BC) 2 - K j F = 0 


(17) 


and a A , b A , a B , bg are complex coeff icients to be determined by 
the boundary conditions. 

f(C) and g(C) have the forms* 
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E c n <i-®c> a . c 0 • i 

n even 

(18) 

£ c (i-B O n , C = 0, C = 1 
n odd 

K C (n-2) " C (n-4) q 9 ) 

n(n-3)B 2 

Hie Boundary Conditions 

At the jet centerline we have already provided for apply- 
ing the condition of zero normal displacement (symmetrical case) 

or zero pressure (anti-symmetrical case). Across each interface 

. ★ 

between two regions the pressure, Ap, and the normal velocity , 
w, must be continuous. (In uniform-flow regions w = 0 z and 
A P = - pf6 t t with coordinates fixed in the fluid. In shear 
layers w is given previously and Ap = p c s.) 

The application of these boundary conditions and the solu- 
tion of the resulting set of simultaneous equations is a tedious 
but straight-forward process. Only the results are given here. 

The Transmission Coefficients and Final Results 

The transmission coefficients may be written 


f(C) = 

g(c) = 

where 

C = 


*More precisely the normal displacement must be continuous, but 
where there are -no discontinuities in streamwise velocity, as 
in this case, normal velocity is also correct. / 
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S" = 


S' = 


-2(1-K M ) 
1 s 


1 r 


_f(0)[3 g -g(0)[ L f c cc 1 ) r 


Den. 


±2(1-1^) '/(l-K^g) 2 -^ [f(0)[ Hg-g(O) [' ] f _ 


Den. 


* I f(C l )[ V g(C l )[ If. 


where in the symmetric case (S’ = S s * , S" = S s ") 

L J g = g ? (C 0 ) cos«r 0 + J[ l+Kj^CMj-Mg)] 2 -K 2 sinc-Q g(? 0 ) 

■ L ]f = f ? (C 0 ) cos^ 0 + /[l+KjCMj-Mg)] 2 -^ sin«r 0 f(( 0 ) 

and in the anti-symmetric case (S’ = S_ ’ , S" = S_") 

r ' 2t 2L 

C J g p sintr 0 - '/Ll + K 1 (Mj-M s )] 2 -K 2 cos<r 0 g(C Q ) 

[]f = f((Co) sin<r 0 - |/ ? [l+K 1 (M j-M s ) j 2 -K 2 cos <r 0 f(? Q ) 

and where 


= 


O' n - 


t(z=z 1 ) = M s a./b 
C(z=z 0 ) = -„(Mj-M s )/b 


«’ z Q //[ 1 + Kj ( M j - M g )] 2 -K 2 " 


Den. = [f c (Ci>L Jg-g/Ci^t Jf 


g * C 


* [a-K 1 H s ) 2 -K 2 | i f <{i)C ig-g( {l )t ] f _ 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 

(24) 


(25) 
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(The sign choice in S’ corresponds to the sign choice in Eq. 
(7).) 

Application of the far-field analysis (using the principle 
of stationary phase) and transformation to retarded coordinates 
gives the critical K values as 


= X 2 /R 

1 1+M s x 2 /R 

= y/ R 

2 1+M s x 2 /R 


(26) 


where x 2 ,y,z (and R = ■ x 2 2 +y 2 +z 2 ) are the retarded coordinates 
fixed in the nozzle*. 

From the far-field potential for a permanent source of 
frequency w (the case so far considered), the far-field poten- 
tial for the transient source sequence described previously 
(Section 2) can be obtained by methods derived in Ref. 4. The 
final result for mean-square pressure in the far field, assuming 
each transient source goes through 2n cycles of frequency <u t is^ 



2A 


2 2 2 7 2 

P "o z 


R 


R 


1+Mx 2 /R 


2 - 0 


GO 

f U/u » 0 ) 2 (S ,2 +S" 2 ) sin 2 (2i7-nw/«» 0 ) d(«/« o ; 


» n 


(u/oi ^) 2 - lj ^ 


(27) 


or, if each transient source is assumed to go through many 
cycles, 


0 2.2 2 2 _2 
~ 2 = 2y A l *0 — 


S ,2 (« 0 ) + S” 2 (u) 0 ) 


(28) 


R 2 R 2 |l+Mx 2 /Rr 

*The nozzle is here assumed to be stationary in the ambient air. 
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In these equations S* = -§-(S s ' +S a * ) , S" = -§-(S s "+S a ") , subscripts 

s and a denoting values found from Eqs. (20) for the symmetric 

and anti-symmetric problems, respectively. The result is the 

mean-square pressure for a single off-center source of strength 
2 

8 ir a. M is the source Mach number (M_). (For a single source 

s 

on the far side of the jet, the mean-square pressure is given 
by Eq. (27) or Eq. (28) with S« = i(S s ’-S a '), S" = i(S s "-S a ") . ) 

Henceforth, for convenience we drop the subscript 0 on the 
generating frequency which will be designated simply w . This 
should cause no confusion. 

Special Case; Source in Central Uniform-Flow Region 

The special case of a source at the center of a shear 
layer jet has been treated in detail in earlier work * . If 
the source is not on the centerline but is in the unifprm-flow 
region inside the shear layers, the analysis is similar. There 
is no need to consider the special 8-region introduced when the 
source is in the shear layer y and the component source solution 
(i.e. the first term of Eq. (6) ) is added to the central region 
potential in Eq. (3) , where now M g = Mj. In the shear layer 
we define 

C = u>’ ( z-Zq ) (29) 

(replacing the previous Eq. (15)), so that in the definitions 
of Eqs. (21) -(25) we have 

«o V ~° 

<r 0 = w'ZqV^L-K 2 


( 30 ) 
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The transmission coefficients reduce, in the symmetric case, to 

3B(1-M s K 1 )V / O^M s K 1 ) 2 -K 5 f ( C x ) L J g -g C C x ) L J f 


S * = +2 COS cr 
s s 


Den. 


S " = 2 cos* 

o S 


3B(l-M s K 1 )[f c (C 1 )L Dg-g^t^L If] 


Den. 


(31) 


where 


o- = u>'z >/l-K 2 ' 
s s 


(32) 


and L J f » [ ]g are defined in Eqs. (21) (with £ ^ = 0). For the 
anti -symmetric transmission coefficients, cos<r g is replaced by 
sin<r g and L j f , L ]g are of course given by Eqs. (22). 


B. DISCUSSION OP COMPUTED RESULTS FOR 
SOURCES OFF THE CENTERLINE 

Typical Examples (for Ms = 0.7) 

Fig. 5 illustrates the effect of moving a source off the 
centerline within the central uniform-velocity region of the 
jet. The greatest mean-square pressure occurs when the source 
is on the ’’near' 1 side (.nearest the side on which the measurement 
is made). However the variation is small, and the average of 
near and far side sources (shown in Fig. 6) is negligibly dif- 
ferent from a source on the centerline. These results are for 
a Strouhal number of 0.2, and would not necessarily hold for 
Strouhal numbers of 1.0 or 2.0. 
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AP 2 R 2 /(2ir 2 A 2 ^ 2 . 2 ) 


Fig. 5 POLAR PLOT OF FAR- FI ELD MEAN-SQUARE PRESSURE 

FOR THREE SOURCE POSITIONS IN THE UN I FORM- VELOCITY 
PART OF THE JET (Two-Dimensional Jet, M ^ = 0.7, 

Srp **• 0.2, z q/ z i = 0.25) 




15 10 5 0 5 


Ap 2 R 2 /(2,r 2 A 2 p 2 «, 2 ) 

Fig. 6 FAR- FIELD MEAN- SQUARE PRESSURE FOR SEVERAL 
SOU ICE POSITIONS IN A SHEAR-LAYER JET 
(Two-Dimensional Jet, Mj = 0.7, S T =' 0.2, 
Zn/ z. = 0.25.) 
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In Fig. 6 the central region of uniform velocity is one- 
quarter of the jet thickness (as in Fig. 5). However the 
sources are placed in the shear layer one-third and two-thirds 
of the way out in the linear velocity profile region. Here the 
far-field mean-square pressures shown are the average for a 
near side and a far side source, whose effects are similar at 
this Strouhal number of 0.2. Fig. 7 shows similar results for 
a source at the mid-point of the shear layer when the central 
uniform-velocity region covers half the jet thickness. 

The most obvious conclusion to be drawn from Figs. 6 and 
7 is that a source of given strength and frequency radiates 
much less noise when traveling at reduced velocity in the shear 
layer than when traveling with the central portion of the jet. 

If the shear layer is to compete with the central region in 
noise production it must contain many more or much stronger 
sources. This is, of course, the case where the potential-core 
region of the jet persists. Further downstream the central 
region of the jet may become important. 

Residual Source and Dipole Radiation in the Present Analysis 

It is also of interest to observe in Fig. 6 the sources one- 
third and two-thirds of the way out in the shear layer, assuming 
them to be 180-degrees out of phase. Such a source-sink pair 
could not radiate at all like a dipole because of the great 
difference in the magnitude of the noise radiated individually. 

The radiation would be primarily that of a source, and this is 
why we stated earlier that elementary reasoning does not appear 
to justify the elimination of source and dipole radiation in 
favor of the less efficient quadrupole radiation in this analysis. 
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Ap^ f?/(2ir 2 A 2 f> 2 '* 2 ') 

Fig. 7 FAR- FIELD MEAN-SQUARE PRESSURE FOR SEVERAL 
SOURCE POSITIONS IN A SHEAR- LAYER JET 
(Two-Dimensional Jet, Mj = 0.7, S^, = 0.2, 

V z i = °* 5 * ) 
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Comparison with Lighthill's Analysis 

The disappearance of source and dipole radiation in 
Ligh thill's analysis is a peculiarity of the highly simplified 
model chosen to represent the jet. In that case the entire air 
mass (jet as well as ambient air) is essentially motionless, 
the acoustic singularities being convected, all at the same 
velocity, through the jet region. This means, for example, 
that any two sources of the same strength and frequency radiate 
energy to the far field in precisely the same amount and pat- 
tern except for phase. Position of the sources in the jet is 
unimportant for this model. Thus cancellation of the source 
strength causes cancellation of the source-like radiation in 
the far field. Similarly, cancellation of the dipole strength 
causes cancellation of the dipole-like radiation in the far 
field, leaving only quadrupole effects. 

In contrast to this model, the present one puts the jet air 
in motion with velocity depending on the chosen velocity profile 
The acoustic singularities travel with the jet air. Then two 
sources, for example, having the same strength and frequency 
generally do not radiate to the far field in the same amount 
and pattern. Consequently, cancellation of the source strength 
does not generally cancel the source-like radiation in the far 
field. 
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5. THE CIRCULAR CYLINDRICAL JET 
AT CONSTANT TEMPERATURE 


A. DEVELOPMENT OF EQUATIONS 


I. THE SHEAR LAYER 

To study the transmission of acoustic disturbances through 

an annular fluid layer in which the velocity (U) varies with r, 

we start with the basic continuity equation and Euler's equation 

7 

(Eqs. 1.2 and 2.3, respectively, in Landau and Lifshitz ) 

|f + div(pV) =0 (33) 

— + (V»grad)V = - i grad p (34) 

dt P 

where p is density, p pressure, and V is the velocity vector. 

In cylindrical coordinates (x,^,,r) we will denote three mutually 
perpendicular velocity components as u,v,w, where u is the x- 
direction velocity (parallel to the axis of the jet), w is the 
radial (r-direction) velocity in a plane normal to the jet axis, 
and v is the tangential (^-direction) velocity in that plane. 

In this notation the continuity equation is written 

+ |U/>u) + I ^-(/> v ) + S~(P Vi ') + ^ = 0 (33a) 

dt dx r r dty r dr r 


and the three components of the Euler equation are 
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du 

at 

av 

at 


at 


• + u— 

+ 

V 

au 

+ 

= - I 2E 


ax 


r 



ar p ax 


. + u^. 

+ 

V 

dv 

+ 

av + vw = . 


ax 


r 

dy 


ar r 

pr di> 

f dvi 

. + u— 

+ 

V 

aw 

+ 

w . ri = - 

l a£ 

ax 


r 

a* 


ar r 

P dr 

u = U(r) 

+ 

Au 

* 

p - p (1+s) 

and 


(34a, b,c) 


O O 

where s is the "condensation" and * = p c (see Lamb , pp. 

o 

476-479). (p Q is the density of the undisturbed fluid in the 
shear layer and p^ is the mean pressure which is constant at 
its ambient value.) Assuming small perturbations and retaining 
only first-order terms, the equations of motion become 


Au t + U(r) Au x + w U r (r) = - c 2 s. 


r t + U(r) v x = - c 2 s^/r 


(35) 


w t + U(r) w x = - c^ s r 


s,. + s v U(r) + Au v + — £ + w_+ll=0 
la a j- r £ 


(36) 


Eliminating Au, v and w in turn gives 
s »t * 3Us ttx * ^txx * U 3 Sxxx 

r 

- C 2 j (v 2 s) t + U(v 2 s). 


Let 


+ 2 V s rx 


(cos) 

s = F(r) cos(m^) sin(kx-wt) 


= 0 


(37) 


( 38 ) 
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where k is wave number in the x-direction and w is the frequency. 
The tangential wave number, m, is integral so that s is single- 
valued. Substituting Eq. (38) into (37) gives 


(1-XJ*MF„ ♦ | ( 2U^k/<n ) + (1-Uk/c)/ (-7+r 0 w' ) | F^ 

+ (1-Uk/u,)! (1-Uk/ W ) 2 - ( k/ <*> 1 ) 2 - m 2 /0 ? +r o u’) 2 jF = 0 (39) 


where tj = «'(r-r 0 ) and u> ’ = u/c . 

For a linear velocity profile U(r) = b(r-r Q ) , and if we 
denote B = bkc/u> 2 = (b/u>)K where K = k/u» * , the differential 
equation of the shear layer is 


(l-Bi 7 )( 7 +r o a.’) 2 F + (i 7 +r 0 «tf’ ) 2B(y+r o «') + (l-Bi;) 


VV 

+ ( 1-Btj) j (ij +r Q «' ) 2 ‘^ 


(1-B„)^ - K z ?- m‘ 


F.= 0 


(40) 


A solution of this equation can be obtained in the form 


F= £-v<i- B ,> n 


(41) 


0 


where, if we denote = 1+Br 


C 0 = Of = ^ = ° 
Go is arbitrary 


(42) 


and 


n(n-3)B 2 u. 2 C = ! B 2 A (n-l) (2n-7) C . 

n i r n “l 


- < (n-2)(n-4)B 2 -K 2 /i 2 -m 2 B 2 > 0^2 

- 2 m K 2 C n _ 3 -(/a 2 - K 2 ) C n _ 4 


'n-6 i 


+ 2 F- C n-5 “ 


(43) 
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(Note that C n is zero if n < 0.) 

To obtain a linearly independent second solution we set 


P 2 = 


F In |1 -Bt, 


w 

E 

0 


A n ( I-B 17 ) 


n 


(44) 


where is the solution already obtained (Eqs. (4l)-(43)) and 

„ 3 2 ,, 2 7 , 2 2 V 

A Q = 3 m B C 3 /(m B -K ^ ) 


A 1 = ° 

A 2 = - 3 M C 3 (B 2 ra 2 +K 2 M 2 )/2(B 2 m 2 -K 2 At 2 ) 


(45) 


l 3 is arbitrary 


and 


n(n-3) /i 2 B 2 A n = 


! - M 2 B 2 (2n-3) C n 


+ llB 


{(4n-9) Cn.j + (n-l)(2n-7) A n _ j ] 
- 2B 2 (n-3) C n _ 2 - { (n-2) (n-4)B 2 

-.K 2 /i 2 -m 2 B 2 |' A n _ 2 

-2K 2 m A n . 3 + (K 2 - M 2 )A n _ 4 


^ 2 ix. A 1 - — . A , 

r n-5 n-o ; 


(46) 


(A n is zero if n< 0.) 

For simplicity we designate the arbitrary constants as a,b 
and denote F^ ~ f , F 2 ~ 8 » and write 


F(y ) = a f ( 17 ) + b g(y) 


(47) 
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with 

f( V ) = £ c n (I" 6 ?)” ' ( 48 ) 

0 

00 

g(^) = f(^) ln|l-B 7 | + A n (1-B^) n (49) 

0 

where = 1, Ag = 1 and C n , A n as above. (b in Eq. (47) 

should not be confused with the slope of the velocity profile 
used earlier.) 

£qs. (48) and (49) express two independent solutions of the 
differential equation (Eq. (40)) as series expansions about the 
singular point -q = l/B. The series converge for those points y 
whose distance from y = l/B is less than the distance from 
y = l/B to the other singular point of the differential equation 
at y = - w 'r 0 . Hence these solutions may not be valid over the 
entire shear layer. In such cases another series solution may 
be obtained by expanding about some other point; we choose the 
outer edge of the shear layer, y = y^ = «'(r^-r o ). Th e con “ 
vergence of this new series is then limited to those values of 
y which are closer to y than the nearer of the two singular 
points of the differential equation. This solution of Eq. (40) 
is obtained in the form 

OO 

P = a n (T 1 ^) 11 (50) 

0 

a^ and a^ are arbitrary and, if we denote (1-By^)/B = * , 


Rj, we obtain 
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r ' n 

nCn-DvRj ^ 2 a n = - . (n-l)R 1 ^(2n-3) + R 1 (4-n) ’a n-1 

+ {(n-2)(n-3)(v-2R 1 ) + (n-2)(3R 1 + l ;) 

+ (Q 0 -i'm 2 ) ; -a n _2 

- < (n-3)(n-5)-(Q 1 +m 2 )|>a n _ 3 + Q 2 a n -4 

^ ✓ 

+ ^3 a n-5 + <Vn-6 “ (51) 


where 

Qq = vR 1 2 (B 2 v 2 -K 2 ) 

Qj = 2R lV 3 B 2 - 3i/ 2 B 2 R 1 2 - 2vR 1 K 2 + K 2 ^ 2 

Q 2 = v 2 B 2 - + 3vB 2 R 1 2 - „K 2 + 2R X K 2 (52) 

rv - a t.2 n -2 2 0 2 n 2 n 2 . „2 

Qj = 6vB Rj - 3» B * B + K 
C ^ = B 2 (3^-2R 1 ) 


(a Q is zero for n < 0 .) 

In this case two independent solutions are obtained by 
letting first one and then the other of the arbitrary constants 
be zero. Thus we may write, as before, 

F(i?) = a f ( 17 ) + b g(y) (53) 

where a,b are arbitrary and 
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In some cases neither series solution (Eqs. (48) and (49) 
nor Eqs. (54) and (55)) is valid over the entire shear layer. 
One series can be applied in one part of the shear layer, 
however, and the other series in another part of the shear 
layer, and there is an overlapping region where both series are 
valid, in such cases the solutions must be matched at some 
point within this overlapping region. This is accomplished by 
establishing conditions of equal pressure and equal radial vel- 
ocity at the matching point, so that the arbitrary constants 
a,b in Eq. (47) are thereby related to those in Eq. (53). 

(The second singular point of the differential equation, 
i.e. that at rj = -w’r D , corresponds to r = 0, the center of the 
jet. As long as the shear layer does not extend in to the jet 
center (a condition which is met by the numerical examples 
studied later) it is not required to obtain a solution by ex- 
pansion about this point. In some cases such a solution might 
be useful to expedite convergence, however.) 

II. THE SOURCE IN A CIRCULAR SHEAR-LAYER JET 
Basic Solution 

In cylindrical coordinates, the wave equation governing 
the propagation of small disturbances in a homogeneous fluid 
is 

I JL(r 22) * i, * 

r dr dr r 2 d if? dx 2 c 2 d t 2 
where 0 is a velocity potential, c is the speed of sound and 


± £& 


( 56 ) 
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t is time. x,f,r are the cylindrical coordinates of a point 
in the homogeneous fluid. 

A solution of Eq. (56) for an acoustical point source, 
located at the origin and at rest with respect to the fluid, is 

f i ( kx- wt ) (1) r — — -v 

0 S = irA R.P. J e H q (rV*' 2 -k 2 ) dk (57) 

-OO 

2 

where A = (source strength)/8ir , source strength being defined 
as maximum volume emitted per unit time, « is the source fre- 
quency, o)' = u»/c , k is wave number in the x-direction, H^^a) 
is the Hankel function of the first kind of order zero and 
argument a, and R.P. denotes "real part of". 

Eq. (57) represents the source as a superposition over 
all wave numbers of cylindrical waves (k «/c ) and exponential 
disturbances (k > w/c). 

If the source is not at the origin, but at some point 
0,^ s » r s* an appropriate wave- expansion formula (see, for ex- 
ample, Magnus and Oberhettinger 9 , p. 21, formula 3b) applied 
to Eq. (57) gives 

CO 

» r .. /* i(kx-u)t) r 

0 = vA 2j « m cos [m(^-^ s )j R.P. J e V 1 **' - k ^ * 

0 —CO 

(1) 2 9" 

• ^ ) dk ; r > r g 


00 

= irA 

. , 0 


•m cos 


CO 

[»(t-* s )] R.P. /'e l(kX " u ’ t) J m (r (6»^-k 2 ) . 


( 1 ) ! 5 2 ' 

• H m (r s V*»* -k^ ) dk ; 


r < r. 


(58) 


\ 
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where «q = 1, « m = 2 for m = 1,2,3, J m is the Bessel func- 
tion of the first kind of order ra, and is the corresponding 

Hankel function. 

Eq. (58) can be derived alternatively by representing the 
off-center source as a 8-function constructed by superimposing 
an infinite Fourier series of ring-sources, each of radius r g 
and having tangentially-varying source strength proportional 
to A cosfm(^-f_)| . The velocity potential of the m th ring- 

L. ^ J 

x. H 

source is given by the corresponding m xn term of the series in 
Eq. (58). (Viewing the source in this manner has practical 
utility later on when computing the far-field pressures.) 

Source in Shear Layer 

To analyze the problem of a source located in the shear 
layer of a cylindrical jet (Fig. 8), it is convenient to first 
assume an annular region (Ar = 28) of uniform velocity, U s , 
within the shear layer. The source is located in this uniform- 
velocity region, and is traveling with it, so that the source 
is at rest relative to the immediately surrounding fluid. In 
a coordinate system moving with the source the jet velocity 
distribution is as sketched in Fig. 9 . This problem can be 
solved, and the thickness of the annular region containing the 
source is then allowed to vanish (8— ►()). 

Since the jet is symmetrical we can expect, on the average, 
to find similar sources at all angles in the jet, but with 
random phase relative to each other. Their mean-square pres- 


sures are additive 


/ 
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r o 

0 


Fig. 8 ST REAM WISE GRGSS-SECTION OF CIRCULAR 

CYLINDRICAL JET WITH SOURCE IN SHEAR-LAYER, 
SHOWING VELOCITY PROFILE AND GEOMETRY 
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Fig. 9 JET GEOMETRY AND VELOCITY PROFILE 
(Coordinates fixed in the source) 
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Coordinate Systems 

The basic coordinate system (x,^,r) is fixed relative to 
the source, with x the streamwise coordinate parallel to the 
jet axis, r the radial coordinate measured from the jet center- 
line in planes normal to the jet- axis, and ^ the angular coor- 
dinate in such planes. Coordinates fixed in the ambient air 
are Xg.^.r and those fixed in the uniform-flow region at the 
center of the jet are Xq,^,T| where 

x~ = x-U t 
2 s 

(59 

*0 = X.(U j -U s )t 

Here U and U. are, respectively, the absolute magnitudes of 
s J 

the source and jet velocities relative to the ambient air, and 
t is time. 


,Tet Temperature and Mach Number 

Most of our previous work has assumed that the jet has the 
same temperature as the ambient air. We now allow the jet tem- 
perature to have a different value from ambient. The jet tem- 
perature, Tj, is assumed constant throughout the jet (r < r^). 
Corresponding to this, the speed of sound in the jet is c^ and 
t,he density is py Outside the jet (r > r^ the ambient air has 
temperature T a , and corresponding speed of sound c a and density 
p & . The temperature ratio is denoted a T : 



( 60 ) 


Although the temperature discontinuity at the outer jet 
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boundary is unrealistic, the analysis will enable a preliminary 
investigation of the effects of heating the jet. 

We define jet and source Mach numbers based on jet tem- 
perature : 
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M s =V c j 

For convenience we define also a Mach number based on the speed 
of sound at ambient conditions so that 


M = U ,/c = '/ST M . 

a; / a ' T j 


M a s * V c a * ■ /i T M s 


The frequency parameter is 

<0 * = ID / c j 

"a* = w / c a = 

A non-dimensional wave number is defined as 


( 62 ) 


( 63 ) 


K = k/«* 

K a = */V = k /^T 


( 64 ) 


The Central Uni form- Flow Region e 

To satisfy the wave equation, provide for outward and 
inward moving waves, arid satisfy the condition of zero radial 
velocity ( d0/di = 0, i.e. no source at r = 0) at the jet 
centerline, 0^ is chosen to be 
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! i »' LKiXo-(U:-U s )t 1-C :tj 

00 (K) = *A«* € m cos|_m«r-# s )] R.P.J % e 


• J m (r«*'/[l+K(Mj-M s )] 2 -K 2 ) 


(65a) 


w - 

0 n = Z / 0 n 00 dK (65b) 

u 0 J u m 

•00 

where ^ is a complex coefficient, = f^* + il^", to be 
determined. 

Eq. (65a) is the component potential for one value of m 
(i.e. a single one of the source rings mentioned previously) 
and one value of wave number K. The total potential is ob- 
tained by integrating over all K and summing over all ra, as in 
Eq. (65b). Since the individual component terms must each sat- 
isfy the differential equation and the boundary conditions, we 
will henceforth write only the component solution for each re- 
gion of the flow field. The total potential is in each case 
the integral (over K) and the sum (over m) of the component 
term, as indicated for the central -region potential in Eq. (65b). 

If the source is located in the central uniform-flow re- 
gion, the source solution from Eq. (58) must be added to Eq. 

(65). For greater generality the following analysis will be 
applied to the case with the source in the shear layer of the 
jet, but the final result for the special case where r g < Tq 
will be given later. 

The Source Region 


In the annular region 8 we have 



- 44 - 


4 <» - 


m 


tA#' « m COS 


[rafy-^)] R.P. j exp [i«' (Kx-Cjt)j 

* { J m ( r s w ’ ) 

(1) / o' _ (2) / 2' n 

+ E* *m <*-’^-K 2 > + G m H. (r«« '/l-K 2 )}J ; 


r_+8 > r > r t 


(66a) 


= tA«’ € m cos[mty-^ s )J R.P.i^exp |_i«’ (Kx-Cjt)| 

. f J (r«* '/l-K 2 ) H^lr «• Vl-K 2 ' ) 

( m m s 


Em H^k.* '£**') + G 2 ' )}j ; 


^ . ( 2 ) 




r g > r > r g -8 (66b) 

The first term in each of these equations is the m,K-component 
of the source potential and the other two terms account for 
outward and inward moving waves reflected off the boundaries 
of the 8-annulus. E m and G m are, respectively, E m ’ + iE m ", 

G m ’ + iG m ", complex reflection coefficients to be determined. 

The first term of each of these equations, integrated over 
all K and summed over all m, gives the source potential which 
at r = r g , </> = <fr s fulfills the necessary conditions for the 
existence of the source. The other two terms, being the same 
in Eq. (66a) and in Eq. (66b), represent a solution which is 
continuous in both pressure and velocity across the 8-annulus. 
Thus the boundary conditions at r = r s are satisfied. 
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The Ambient Air 

Outside the jet (Region 2, the ambient air, in Fig. 9) 
we have 


0_ (K) = tAw* « cosTmC^-^ )j R. P. 
2 m di L s — * 


J m 


H^r«> vTcT7c a ) 2 (l-KM s )2-K2') 


exp^ioi' | K(x 2 +U g t)-Cjt j] 

(67) 


where S m is the complex transmission coefficient, S m = S m ’ + iS, 


m m 


m 


Note that = «/cj and K = k/w* , as before, so that the c^ mul- 
tiplying t in the exponential term is due to this definition of 
«' in terms of jet temperature rather than ambient. The temper- 
ature difference effect appears in the argument of the Hankel 
function. That is, Eq. (67) is a solution of the local wave 
equation (Eq. ( 56 ) t with x = x 2 and c = c a ). 


When KM > 1 and (c-/c ) 2 (1-KM ) 2 > K 2 , nf 1 } ) in Eq. ,(67) 
5 J £1 s m 


( 2 ) 


must be replaced by H ra ( ), the Hankel function of the second 
kind. This is, in cylindrical coordinates, the ’'reversed wave” 
case discussed in previous reports 2 . 


The Shear Layer Regions 

So far we have considered only the regions of uniform flow. 
The shear layer regions must be treated differently, using s, 
the condensation, and w, the radial velocity, instead of the 
velocity potential 0. The development of the shear layer equa- 
tion for a circular cylindrical jet with linearly varying velo- 
city profile was given earlier. As shown then, the conventional 
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wave equation is replaced by a more complicated form which, 
upon the assumption of simple harmonic forms in the x and 
directions, yields an ordinary differential equation (Eq. (40)) 
in the r (or 17 ) direction. Two independent solutions ( f ( 17 ) and 
g(y)) of the shear layer equation are found by series expansion 
and have the forms shown previously. These are combined to give 
the shear-layer region solutions used below. 

In Region IB (see Fig. 9 ) we may write 


m 


w 


m 


cosj m(^-^ s )j 
cos[m(*-^ s )j 


R.P. 

R.P. 


W e 


iu»« (Kx-Cjt) | 


l. m 


1C; - 


J_Q 


i«’ (Kx-Cjt) 


•B n 


1-Bij b m-i/ 


( 68 ) 

(69) 


where B = (b/w)K and b = du/dr. U is the local shear velocity 
in the layer. 17 g is a coordinate in the r-direction defined by 

V B = «»[r - (r s + 8 )] (70) 


Qg is a complex function of i 7 g and is a combination of the 

solutions f m (i 7 ) , g m ( t ?) of the ordinary differential equation 

which applies in the shear layer. Q is the derivative of 

°m»7 

Q D with respect to . 

B m B 

Similarly, in Region 1A (Fig. 9 ) we write 


s m = c ° s i m )j R.P. 


e 


i»' (Kx-Cjt)| 


(71) 


w„ 


= cosj ra (*-*_) I R.P. 


1C; - 


Qa e 
l-Bi7 A ^ i v 


i «' (Kx-cjt) 


( 72 ) 


where 17 is a coordinate in the r-direction defined by 
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^ A = “'fr ' (r S - 8) i 
Now if 8 —0 

^A = ^B = ^ 


i> = «' (r-r s ) 
We may write 



l m 


f n (,,) * b A 


m 





B 


m 




+ b 


B 


m 




(73) 


(74) 


(75) 


where, as stated already, f m ( 17 ) , g m (iy ) are solutions of the 

shear-layer equation. a a , b. , a n , F n are complex coef- 

A m *% °m 

ficients to be determined when all the boundary conditions 
are satisfied. 

The Boundary Conditions 

At the jet centerline the condition of zero radial velocity 
has already been provided for and, similarly, at the radius of 
the source location, r = r g , the solution was constructed to 
provide continuity except for the required existence of the 
source at r = r s , $ = ^ s . Across each interface between any 
two of the flow regions the pressure (Ap) and the radial velo- 
city (w) must be continuous. 

In the uniform-flow regions w = and Ap _=_ .«*_ />.0. if the 

uniform-flow region is within the jet or Ap = - ^ a 0 a outside the 

jet, with coordinates fixed in the local fluid in each case. In 

2 

the shear-layer regions w is given previously and Ap = ^j c j s * 
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Application of the pressure and velocity boundary con- 
ditions (at the four boundaries r Q , r s -8, r g + 8, r i) results 
in a set of eight simultaneous equations to be solved for the 
eight unknown coefficients (a^, 5^, a^, b^, R m , S m , E m , G m ) . 

Solution is tedious but straight-forward. 


The Transmission Coefficient 

For convenience, the following notation is used: 


V 0 = n(r 0 ) = «'(r 0 -r s ) 

*7 1 = = w’Crj^-Tg) 


(76) 


o - 0 = «*r 0 /./[ l+K(Mj-M s )J 2 -K 2 

o-g = «'r s 'Jl-K? (77) 

< r 1 = “' , r 1<# /(c j /c a )2(i-KM s )2-K2' 


*/®' = '/[l + K(Mj-M s )J 2 -K 2 ' 

= , /(Cj/c a ) 2 (l-KM s )^-K 2 


(78) 


Then the complex transmission coefficient is 

2i (Fj//> a > f m^l)gm,^l) - gm^l) f m CJ N 

g = 

(1-KM ) f (0)g m (0) - g(0)f (0) L ] n 

s s m m m ^ D 


where 


t] » = V 0)I Jg m0 - **«»- h 


mO 


CD n = - 3 f C -I _ 

u r m0 g ml 


i j 

g m0 ^ml 


( 80 ) 
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and 


”f m0 = W f n , ( ->0> - ^ V<'o> f m<’o> 


rj g n0 ' J » ( '0 >B «, ( »0 > • ''© W^V 


( 81 ) 


C] f m i = V ( 'l> f .» ,'n> - ^ > <' 1 >V’> 1 > 


( 1 ) 


■x ( 1 ) 


[ ] 


(1) r^s (1) 

gjnl = "»> - <8j/a a )-/© H, •Ca 1 )g B1 ( n ) 


(82) 


The prime following the Bessel or Hankel function indicates 
the derivative (e.g. Jm’Cos) = 

In the reversed wave case (K >• l/M s and <r^ real), H^^o^), 
H^ 1), («r 1 ) are replaced by H^ 2 ^), H^'fo^) in Eq. (82). 


Special Case: r g < r Q 

If the source is located in the central uniform-flow region 
rather than in the shear layer, the preceding analysis is slight- 
ly simplified. In this case we define 

V = u’(r-r 0 ) (83) 

Then »7 x = ^(r^rQ), r, Q = 0; also o- Q = w'^VT-K^. The 
transmission coefficient then reduces to 



2i( />•//> )J (<r ) 
J a - m s 

(v<» ' r Q ) (l-KM g ) 




m ± m ^ 1 


r 

L 



( 84 ) 


(Note that M_ = M ; in this case of course.) 

O J 
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Additional Interfaces and "Matched" Solutions 

Note that we have considered the shear layer to be linear, 
with uniform slope throughout. It is apparent that a more com- 
plex multi-slope profile, as sketched in Fig. 10a, could be 
studied by the same methods if conditions of continuous pressure 
and normal velocity are applied at each interface between linear 
segments of the shear-layer velocity profile (e.g. at r = a, 
r = b, in Fig. 10a). Each interface thus adds two equations to 
the set of simultaneous equations which must be solved for the 
transmission (and other) coefficients, and the form of the trans- 
mission coefficient is then somewhat more complicated than that 
shown in Eqs. (79) and (84). 

Essentially the same method is employed when no single 
series-expansion of the shear layer solutions f(y), g(jj) exists 
which is convergent over the entire shear layer (see paragraph 
following Eq. (55)). In this case the shear layer is divided 
by an imaginary interface at r = r fi , though the slope of the 
shear layer may be the same on either side of r g (Fig. 10b). 
r g is chosen anywhere in the region where both pairs of series- 
expansion formulas are valid; then one pair (Eqs. (48) and (49)) 
is useful in one part of the shear layer (tq 4 r 4 r e ) and the 
other pair (Eqs. (54) and (55)) is useful in r e 4 r 4 r^. Con- 
ditions of continuous pressure and continuous radial velocity at 
the imaginary interface r = r g add two equations to the set of 
simultaneous equations containing the unknown amplitude coef- 
ficients. Two forms of the transmission coefficient are ob- 
tained (for the single-slope velocity profile case) depending 
on whether the source radius r s is greater or less than r e . 

The solution is omitted here for brevity. 
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Par- Field Result and Extension to Transient Source Case 
The velocity potential outside the jet is 

ao 

K'tf 0 2 (K) dK (85) 

0 - OO m 

(67), contains the transmission coefficient S m , 
the determination of which has just been described. In the 
far field, the integral of Eq. (85) can be evaluated by the 
method of stationary phase. This analysis has previously been 
given in some detail for the two-dimensional jet case 2 * 2 , and 
only the results for the circular cylindrical jet will be 
summarized here. 

For '/ap M s (- M ag ) < 1 there is only one point of sta- 
tionary phase, located at K = K c : 

K c = V«T (x/R*-M s /a T )/( 1-M s 2 a x ) (86) 

where 

R* = x2+(i-M s 2 o T )r2' = '/x 2 +(l-M as 2 ) r 2 ^ (87) 

In the far field, then, 


0 2 —-2UA/R*) « m cos[m(*-* s )j 
0 


R.P. 


* ^m^c^ e 


i ( w 1 y c - «t -m *■/ 2) j 


where 


1 r c - <o ’ V j R* - M s v a jXj /( 1 - M s 2 o t ) 


( 88 ) 


(89) 


Using Ap 2 = /j & 2 02^^ (with x = constant in taking the t deri 
vative) gives the mean-square pressure in the far field: 



-53- 


~ 2*2a 2 /> a 2 w 2 


Ap 


R* 2 (1-M s 2 o t ) 2 


M s x 

R* 


12 


00 00 


EE * m *n cos[m(*-* s )] 
0 0 


• cos[n(<j/-^ s )j 


( Sju ' S n ’ +S m "S n ") [cos (^) cos(^) +sin(5*)sinX£pj 


+ (s m' s n”- s m' ,s n , >j cos(^)sin(^)-sin(^)cos(^)j 


IIttn 


rmir> 


r mrx! 


1(90) 


(n is used here and in Eq. (93) as a dummy index in the double 
summation.) 

The analysis so far has been applied to the "permanent" 
source, i.e. an acoustical source drifting downstream with the 
local fluid for all time. Por the sequence of transient sources 
described earlier Ap2 is obtained 2 *'* by formally converting Eq. 
(90) to retarded coordinates and multiplying by (l+M^x^R). 

The retarded coordinates are X 2 »r,^,; R where 

R = \/x2 2+ r 2 

Hie original coordinates fixed in the source are x,r,^;R* 
where R* is defined by Eq. (87). The equations for coordinate 
transformation, from Refs. 3 or 4, are 


R* = R + M a x 2 
i s “ 

x = x 0 + M a R 
z a s 


(91) 


Changing to retarded coordinates the critical K values become 

K c = y ; «f (X2/R)/(1 + Af m s x 2/ r ) (92) 

or, in the notation of Eq. (64), 


K a c = (x 2 /R)/(l +M a s x 2 / R) 


(92a) 
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XpZ for the transient source sequence is 




0 2,2 2 2 
2 *• A p. 


00 CO 


R 2 jl + M as X2/Rj J 0 0 


~3 ?? Vn cos[m(^-^ s )] cos~n(*-*.)j 


( S ra * S n 1 +S ra f?S n n) [_ cos (^) cos (^) +s in (^) sin (4^) 


+ (S ra ’S n "-S m "S n , )rcos(^)sin(^)-sin(5J)cos(|l)! | 

L _i 

(93) 

Note that each transient source is assumed to go through many 
cycles. The transmission coefficients are evaluated at K_ 

c 

• I ! 

(Eq. (92)). When the absolute magnitude of the jl+M as X 2 /Rj 

factor is used, as indicated, Eq. (93) is applicable for M ag % 1. 

Eq. (93) is the far-field mean-square pressure for one 
source (i.e. one transient source sequence) located at </» = $ s , 
r = r g . Since the jet is symmetrical we can expect, on the 
average, to find a similar source located at any around the 
jet, each with random phase relative to all of. the others. 

Their mean-square pressures are additive. The average mean- 
square pressure due to a ring of random-phase sources is 

.2ir 

= JL / df, (94) 


or 


Ap 2 r2 


1+M as x 2 /R J 0 


£ •m<S m ’ 2 +V ,2 > 


(95) 
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Note that each value of m in Eq. (95) represents the 
average mean-square pressure in the far field due to a tran- 
sient ring-source having source strength proportional to 
cosj m(<i»- ^» s )j . Since each of the ring sources with m > 0 has 
zero net source strength, it might be expected that radiation 
from positive sections of the ring would tend more and more to 
cancel the radiation from negative sections of the ring, so the 
net acoustic radiation into the far field would be small for 
large m. Computations for subsonic jet examples indicate this 
to be the case (though it may not be true for supersonic jets). 
tVhen Mj <r 1 the terms of Eq. (95) tend to drop rather quickly 
in magnitude as m increases from zero, so that in practice only 
a few terms in the summation usually are required, and computing 
time is thereby benefited. 

B. RESULTS FOR 

JET TEMPERATURE = AMBIENT TEMPERATURE 

This section contains a comparison of far-field mean-square 
pressures produced by sources in circular cylindrical jets and 
in two-dimensional jets, at subsonic velocities. Comparisons 
are made for the following conditions: 

1. Source at the center of a high Mach number subsonic 
jet. 

2. Source in the shear layer of a high Mach number 
subsonic jet. 

. Source in an extended constant-velocity region 
centrally located in a high Mach number subsonic 
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jet. 

4. Source at the center of a low Mach number jet. 

For the two-dimensional jet, pressures were evaluated only 
in streamwise planes containing the source and normal to the 
plane of the jet. A simulation of circular jet effects from 
such data would at least require evaluation of pressures in 
planes other than normal to the two-dimensional jet plane. Some 
of the differences between circular and two-dimensional jet ef- 
fects pointed out in the following pages can in a crude sense be 
attributed to the neglect of non-normal planes in the two- 
dimensional calculations. 

1. Source at the Center of the Jet; M j = 0.7 

In Fig. 11 the far-field mean-square pressures are shown 
for a source at the center of a circular cylindrical jet. The 
maximum Mach number in the jet is 0.7 which is held constant 
throughout a small region near the jet center. The radius of 
this region is 0.1 times the jet radius, and the jet velocity 
varies linearly outside this region to a value of zero at the 
edge of the jet. The numerical calculations are simplified by 
having a central constant-velocity region, and the cross- 
sectional area of the central region is in this case only one- 
percent of the total cross-sectional area of the jet. Strouhal 
numbers (S-jO of 0.1, 0.2, 0.5 and 1.0 are shown, where 

<,_(*) (jet diameter) 

- - • 

L 2n (maximum jet velocity) 
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Fig. 12 FAR- FIELD MEAN-SQUARE PRESSURE DUE TO A 

SOURCE AT THE CENTER OF A TWO-DIMENSIONAL JET. 
Mj = 0.7; z Q = 0. 
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For comparison Fig. 12 shows far-field mean-square pres- 
sures for a source at the center of a two-dimensional jet. 

Again the maximum Mach number of the jet is 0.7. The jet velo- 
city varies linearly from the jet center to the edge, although 
for strict comparison one should perhaps maintain a central 
constant-velocity region with thickness one-percent of the jet 
thickness. 

Comparison of Figs. 11 and 12 shows that for the circular 
cylindrical jet the pressure lobes are wider and are rotated in 
closer to the downstream axis. For a Strouhal number of 0.2 
the mean-square pressures are of similar magnitude. However the 
variation of mean-square pressure with Strouhal number is con- 
siderably greater for the circular jet. 

On Figs. 11 and 12 (and most following figures) the radial 
line with e on one side and w on the other indicates the boundary 
between the region in which the disturbances leaving the source 
might be termed "pseudo -sound" (e) and are exponentially de- 
caying in z (or r) and the region in which disturbances leaving 
the source are "true sound", or true waves (w). It is evident 
that most of the far-field noise for the examples shown in Figs, 

11 and 12 is produced by disturbances originating as pseudo- 
sound. 

All of the pressure lobes obtained in this study are 
closed curves, i.e. Ap 2 -*0 as ©-»- 0 and as © -*•* . The plotted 

curves have sometimes not been completed in the vicinity of 
© = ir (as, for example, the S T = 0.1 curve in Fig. 11). This 
is because the computer output was obtained only for certain 
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specified ©-values which in some cases were not closely enough 
spaced to allow good fairing of the Ap^-curves near © = v. 

This is a small region, however, and we have simply left that 
part of the curves unplotted, since it did not seem worthwhile 
to re-run the program with more ©-values. 

2. Source in the Shear Layer 

Fig. 13 shows, far-field mean-square pressures for a source 
half-way out in the shear layer of a circular cylindrical jet. 
The source therefore travels at half the maximum jet velocity . 
and its Mach number is 0.35, while the small constant-velocity 
region at the center of the jet has a Mach number of 0.7. 
Comparing this with Fig. 14 , a corresponding case for the two- 
dimensional jet, we see again that the pressure lobes are 
broadened and extend in closer to the downstream axis. For the 
lower Strouhal numbers the variation of mean-square pressure 
with Strouhal number is increased as before. For = 1 and 
.S,p = 2 the two-dimensional jet shows a reversal in trend with 
the peaks moving out from the downstream axis. This reversal 
does not appear for the circular jet, though it might occur at 
still higher Sj. values. Such effects might appear as the 
wavelength of the sound is reduced to the order of the jet 
diameter. For the smaller Strouhal numbers, the far-field 
noise again mainly originates at the source as pseudo-sound, 
but this is less true as the Strouhal number increases. 



3 


2 


Ap 2 R 2 , 


FAR- FIELD MEAN-SQUARE PRESSUR] 
OUT IN THE SHEAR LAYER OF I 
Mj = 0.7; M s = 0, 
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3. Source in an Extended Constant-Velocity Region 
Centrally Located in a High-Velocity Jet 

For the example chosen the constant-velocity region in the 
circular jet extends from the center to one-half the total jet 
radius, and in the two-dimensional jet from the center to one- 
quarter of the jet semi-thickness. In each case one-quarter of 
the cross-sectional area of the jet is in the constant-velocity 
region. 

In Figs. 15 and 16 mean-square pressures (in the far field) 
are shown for sources at the center and at the edges of the 
constant-velocity region. The maximum jet Mach number is 0.7 
and the Strouhal number is 0.2 in all cases. For the circular 
jet (Fig. 15) a source at the center produces slightly lower 
pressures than sources at the edge of the constant-velocity 
region. This occurs in the angular range where source distur- 
bances decay exponentially in the central region (angles (©) 
greater than 126-degrees), and it is probably this exponential 
decay which weakens the effect of the source at the center. 

The two-dimensional jet, Fig. 16, shows a similar result, the 
effectiveness of the sources progressively decreasing as their 
disturbances have to traverse increasing thicknesses of the 
central region. 

The pressure lobes for the circular jet are somewhat 
broader and rotated in slightly towards the downstream axis by 
comparison with the two-dimensional jet. At this one Strouhal 
number of 0.2 the pressures are appreciably lower for the 
circular jet. 




Fig. 15 MEAN- SQUARE PRESSURE IN THE FAR FIELD DUE TO A SOURCE 
IN A CIRCULAR JET, FOR TWO SOURCE LOCATIONS: 
r g = 0 and r g = 0.5. 

Mj = 0.7; S T = 0.2; r Q = 0.5. 
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4. Source at the Center of a Low Mach Number Jet 

Figs. 17 and 18 show far-field mean-square pressures for 
sources at the center of a Mach 0.3 jet. In the two-dimensional 
jet (Fig. 18) the velocity falls off linearly from the center to 
the outer edge. In the circular jet (Fig. 17) the velocity is 
constant out to 0.1 times the jet radius, and then falls off 
linearly to the outer edge. (Only one-percent of the jet cross- 
sectional area is in this central region.) 

As in the higher Mach number case, the principle effects 
of changing from the two-dimensional to the circular jet are, 
first, that the pressure lobes broaden and extend in closer to 
the downstream axis and, second, that the pressure variation 
with Strouhal number is accentuated. 

Note again the radial line marking the boundary between 
the region where disturbances leaving the source are pseudo- 
sound (e) and the region where disturbances leaving the source 
are true waves (w). At a Mach number of zero (a trivial case 
not shown here, with Ap^ r2/ (2-n^A^p^ui^) = 1) all far-field 
noise is produced by disturbances which leave the source as true 
sound. It is therefore remarkable that at a Mach number of only 
0.3 most of the far-field noise is produced_by disturbances 
originating as pseudo-sound. Since pseudo-sound corresponds to 
the 90-degree out-of-phase components of pressure and velocity 
it transmits energy only by interference with reflected waves, 
in the transmission through the jet these exponentially decaying 
disturbances become true waves, and so can transmit energy to 
the far field without assistance. 
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C. RESULTS FOR 

JET TEMPERATURE ABOVE AMBIENT TEMPERATURE 

As an example for numerical computations, we have chosen 
a circular cylindrical jet with uniform velocity throughout 
a core covering the inner half of the total radius. The velo- 
city decreases linearly over the outer half radius, reaching 
2ero at the outer edge of the jet. All comparisons are made 
at fixed jet velocity , so that the jet Mach number based on 
the ambient speed of sound is 0.7 (M_: = 0.7). (Qualitatively 
similar results, not shown here, are obtained at a Mach number 
of 0.3.) The ratio of jet temperature (uniform across the jet) 
to ambient temperature is designated by a^, and values of 1, 2, 

4 and 6 are considered. 

Figs. 19 - 22 show results for Strouhal numbers of 0.2, 

0.1, 0.01, and 0.001, respectively. The very low Strouhal 
numbers are included to show that in the limit of zero frequency 
(or zero jet diameter) the temperature effect disappears. The 
curves then merge and the value of Ap2R2/ ( 2*2^2 ,, 2 W 2) becomes 
jl+M cosej - ^, which is the value for a modified moving source^. 
At higher Strouhal numbers (up to and including 0.2, at. least) 
increasing the jet temperature reduces the magnitude of the 
pressure lobes and rotates them away from the downstream jet 
axis. This trend may alter at still higher Strouhal numbers, 
which have not yet been investigated. 

It must be remembered that the source strength here is 
8 tt^A, and represents the maximum volume introduced by the 
pulsating source per unit time. Only if this source strength 
is unaffected by heating the jet (with velocity profile fixed) 
will the reduction in mean-square pressure with increasing 
jet .temperature be observed. 
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Fig. 19 PAR-FIELD MEAN-SQUARE PRESSURE DUE TO SOURCE 
AT CENTER OF UNIFORMLY HEATED CIRCULAR JET. 

Ma j = U j/ C a = ® #7; = °* 2; V 0,5 ‘ 
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Fig. 22 M a j = Uj/c a = 0.7; S T = 0.001; r Q = 0.5. 

FAR- FIELD MEAN-SQUARE PRESSURE DUE TO SOURCE 
AT CENTER OF UNIFORMLY HEATED CIRCULAR JET. 
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6. PRELIMINARY ANALYSIS FOR 
LINEARLY VARYING TEMPERATURE IN JETS 

Introduction 

In previous analyses we have considered mean* velocity 
variations across the jet, but the mean temperature of the jet 
has been held constant (either at the ambient value or higher). 
Here we permit temperature to vary across the jet also, and car- 
ry the development as far as the ordinary differential equation 
for the shear layer in both two-dimensional and circular jets. 

For simplicity we show the development only for the two- 
dimensional jet. The ordinary differential equation for pres- 
sure as a function of z has one singular point if the tempera- 
ture is constant. This singular point may be in the region of 
interest (the shear layer) or outside it. In the former case 
the expansion must be made about this point and is everywhere 
convergent. Linearly varying temperature introduces a second 
singular point where T = 0. This point never falls in the re- 
gion of interest, but it does limit the radius of convergence 
for the expansion about any other point. In some cases two 
matched expansions are then required, which complicates the 
computation. 

In a circular jet a third singular point appears (at the 
center of the jet) and this, imposing further limits on con- 
vergence, might require more than two matched expansions in 
some cases. Thus the programming of the calculations would be 
more complicated for the circular jet. 

*mean refers to a time average eliminating turbulent velocities. 
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Development of Equations for Two-Dimensional Jet 

Let x,y,z be coordinates measured longitudinally, laterally 
and normal to the "plane" of the two-dimensional jet, and u,v,w 
be the corresponding velocities. Then, indicating perturbation 
quantities (which are functions of x,y,z and t) by Au, etc., we 
have 

u = U( z) + Au 
v = Av 
w = Aw 

T = T q Cz) + AT 
P ~ P 0 (z) + A p 
p - p Q + Ap 

where 

p 0 = constant = Pambient = ^ 0 (z)R T 0 (z) 
and where R is the gas constant, assumed fixed. U(z) gives the 
mean velocity profile, T Q (z) the mean temperature profile and 
p 0 (z) the mean density profile. The mean pressure, p Q , is 
constant throughout the jet and ambient fluid, and the mean 
temperature and density profiles are simply related through the 
equation of state. The dynamical equations can now be written 
as 8 

Ap 

Au + + U(z) Au„ + Aw U_(z) + — t— r* = 0 
T x z Po (z) 

Ap v 

Av^. + U(z) Av x + = 0 

'o (2 > 

AP Z 

AWj. + U(z) Aw x 4- = 0 

/»o (z > 


( 96 ) 
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and the equation of continuity becomes 

A*> t + U(z) A/» x + ^> 0 (z) jAu x + Av y + Aw z J + p Q Aw = 0 (97) 

1 — * — z 


Here we have four equations and five unknowns: Au, Av, Aw, 

Ap and A p. However, if we assume that the entropy of. each fluid 
element is unchanged when an acoustic wave passes through (e.g. 
see Blokhintsev 1 ® , NACA TO 1399, p. 21) then, following a fluid 
particle, 


d(Ap) 

d(A/>) 


k(z) r ^ 0 ( z )] r " 1 = 


y RT 0 (z) 


c 2 (z) 


(98) 


where k(z) is the adiabatic constant for a given z and y is the 
ratio of specific heats. (Note that entropy is not assumed 
constant throughout the jet nor even at one point in the jet.) 
c(z) is the local speed of sound. 

Prom Eq. (98) 

d A p _ dAp 1 (99 

dt dt c^(z) 


Also, the substantial derivatives, dA^/dt and dAp/dt, are 


= A/> t + A , o z (z) Aw + U(z) A f>x 

dAp 

= Ap t + U(z) Ap x . .i; • 


( 100 ) 


Eliminating dAp/dt and dAp/dt between Eqs. (99) and (100): 

A p t * U(z) A fix = ;^P t + U(z) APxj ”/o z (z) Aw (101) 

This relation between Ap and A p gives us the required added 
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equation to make the previous set determinate. If, for con- 
venience, we express Ap as 

Ap = /> 0 (z) c 2 (z) <r (102) 

(noting that /> Q (z) c 2 Cz) is a constant), then the four equations 
become 

Au t + U(z) Au x + U z (z) Aw + C 2 (z) <r x = 0 

Av, + IKz) Av + c 2 (z) a = 0 
t x y 

(103) 

Aw^. + U(z) Aw x + c 2 (z) <T Z ~ 0 

+ U(z) <r x + Au x + AVy + Aw z = 0 


Eliminating Au, Av and Aw gives 
*ttt + 3 U(z) °ttx + 3 u2tz) ' 

- C 2 (z)| (V 2 <t) + U(z) (V 2 <r) 

I x 


txx 


* u 3 <z) ' XXX 

+ 2 U z (z) c 2 (z) 


<r 


zx 


- [c 2 (z)] z hr 2t + u(z) oj = 0 (104) 

For [c 2 (z)J z = 0 (i.e. a constant temperature jet) <r be- 
comes s, the condensation, and Eq. (104) becomes Eq. 3 of Ref. 4 
or Eq. 5 of Ref. 1. 

For U(z) = 0, the case of a static heated layer, we get 


V 


2 

<7 = 


c 2 (z) 


tt 


c 2 (z)] z 



(105) 
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Re turning to Eq. (104), let 

cr(x,y,z,t) = F(z) cos(k 2 y) siiHkjX-wt) 


Substituting this in Eq. (104) gives an ordinary differential 
equation for F as a function of z. 


F 


zz 


+ 


2U z ki ^ (c 2 ) z 

+ P 

— — ■ 

<„-Uki) 2 ,_ 2 

U-1%) c 2 

“ r 

1 

[ 

O 

CO 

1 

► 

* 

I 


0 


(106) 


The simplest temperature (or c 2 ) variation with z is 
linear. In that case the equation has two singular points: 
one, where (w-Uk^) = 0, may fall within the region of interest 
(the linear shear layer); the other, where c 2 = 0, never lies 
in the region of interest. If F is expanded as a power series 
in z, one expansion must be made about each singular point in 
the region of interest. Other singular points limit the radius 
of convergence of such expansions and may require the use of 
matched expansions about two different points. Such procedures 
have already been described in connection with the circular jet 
at constant temperature, which has two singular points. 


The Circular Cylindrical Jet 

For the circular cylindrical jet the axial, angular and 
radial coordinates are x, f and r. The corresponding pertur- 
bation velocities are u,v,w. 

The partial differential equation for o-is 
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- C 2 KV 2 <r) t + U (V 2 <r) x ; + 2 IJ r C 2 <r rx 

- (c 2 ) i<r +U«r j = 0 (107) 

r' u rt irx_j 

p o 

where U = U(r), c = c (r) and the perturbation pressure, Ap, is 
Ap = /> 0 (r) c 2 (r) <r (/» Q (r) c 2 (r) is a constant). 

If we let 

<r(x,y/,r,t) = F(r) cos(m^) sin(kx-wt) 


then F(r) is defined .by 


rr 


+ F. 


2U r k 

(w-Uk) 


1 d(rc 2 ) | 
rc 2 dr i 



(108) 


Note that k here corresponds to k^ in the two-dimensional jet. 

For temperature (or c 2 ) varying linearly with r, singular 
points occur for ( W -Uk) =0, r = 0 and c 2 = 0. Obtaining solu- 
tions by series expansions is therefore more complicated than 
in the two-dimensional case which has only two singular points 
instead of three. 
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7. DISCUSSION OF APPLICABILITY 

As stated earlier this work is not intended to provide 
a comprehensive theory of the noise produced by turbulence in 
jets. Instead, we have concentrated on one aspect of the prob- 
lem with the purpose of providing facts (not intuitively evi- 
dent) which can be used in formulating a more comprehensive 
theory. Utilization of this material in a more general theory 
requires consideration of a number of factors. Some of these 
(not necessarily in order of importance) are as follows: 

(1) What singularity (source, dipole, quadrupole or other 
type) should be considered as the noise generator? 

(2) Since we consider that turbulence is the only true 
generator, we are interested in experimental methods capable of 
distinguishing between turbulent velocity fluctuations (asso- 
ciated with vorticity and non-linear terms) and the linear 
fluctuations associated with transmission in our analysis. 

(3) Frequencies recorded by an instrument at rest in the 

ambient fluid are not the same as the generating frequencies 
(because of the Doppler effect). The frequency correction is 
easily made, but the plotting of mean-square pressure curves 
for fixed observed frequencies requires the use of an energy 
spectrum, since a given generating frequency produces differ- 
ent observed frequencies at different angular positions in the 
far “field."' ‘ ' " ' " " _ 

(4) Various generating regions of the jet must be con- 


sidered 
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(5) The limitations imposed by considering the jet to be 
infinitely long should be studied, especially at low frequencies. 

The final factor ((5)) is now considered in some detail. 
There are really two parts to this question. The first part 
might be stated as follows: Assuming coordinates fixed in the 

source, how much of the jet upstream and downstream of the 
source is really essential to the analysis? The second part of 
the question relates to the transient source analysis in re- 
tarded coordinates: How far does a given transient source tra- 

vel relative to the nozzle, and can it be assumed to pass through 
many cycles in this lifetime? 

The first part of the question relates to refraction ef- 
fects, which do not appear explicitly in Lighthill’s work. 
However, the second part of the question relates to convection 
effects and so may apply also to Lighthill's analysis. 

Part 1 (coordinates fixed in the source) 

When the generating frequency is low the wavelength of the 
disturbance is large, and may be many times the jet thickness. 

Our analysis deals primarily with wavelengths in the streamwise 
direction and, since we assume an infinitely long jet, many 
cycles may appear in the jet for any given wavelength. On the 
other hand, in a real jet of finite length, perhaps only a 
fraction of a cycle may appear. This raises a question as to 
the applicability of our analysis at low frequency. 

However, in certain cases we have examined the results 

obtained for infinitely long jets in the limit of zero fre- 

5 

quency. For example a source drifting with the fluid in an 
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inf initely long jet becomes a "modified moving source"* in this 
limit. Such a result is also what one would anticipate for a 
source in a finite length jet, since the modified moving source 
merely requires a small neighborhood of fluid traveling with 
it. A second example^ is a source in a jet, but at rest rela- 
tive to the external air. This source is necessarily moving 
relative to the jet fluid, and is considered a modified moving 
source. Analyzed by infinite jet methods in the limit of zero 
frequency the result is a simple source radiating uniformly in 
all directions. This is again quite consistent with what one 
would expect for the finite jet in the zero frequency limit. 

Another intuitive argument is obtained by approaching the 
finite jet problem (qualitatively) by modifying the boundary 
conditions for the infinitely long jet. Consider a jet of 
circular cross-section and, with the source as a center, con- 
struct a spherical surface which passes through the jet about 
one and one-half diameters upstream and downstream from the 
source. The portion of the sphere within the jet intercepts 
a solid angle which is five percent of the total for the sphere. 
It is within this limited angular range that the boundary con- 
ditions would have to be modified for the finite jet if uniform 
conditions are assumed over the three-diameter length. This 
suggests, though it certainly does not prove, that the radia- 
tion field outside the jet might not be seriously altered by 
the restriction of jet length. 

It should be mentioned here that there 'is a region in 

which infinite-jet theory is definitely incorrect. This is a 

*The modified moving source is defined and described in Ref. 5 
and compared with the conventional moving source. The existence 
of other types of moving sources also is noted. 
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relatively small angular region in the far field near the jet 
axis. Here infinite-jet theory predicts mean-square pressures 
approaching zero at the axis, which would not be true for a 
real jet of finite length. Since the region is small and the 
pressures are ordinarily small we are not usually much con- 
cerned with this difficulty. 

The preceding discussion may serve to indicate that 
intuitive arguments sometimes conflict and must be handled 
with caution. The range of conditions under which infinite- jet 
analysis is useful must be considered unknown at present. 

Part 2 (retarded coordinates and transient sources) 

The mean-square pressure for a sequence of transient 
sources contains an integral of the sum of the squared trans- 
mission coefficients over the entire frequency range (Eq. (27)). 
There is, however, a weighting factor in the integrand which, 
if each transient source goes through many cycles during its 
existence, becomes a delta-function. In this limit the mean- 
square pressure formula is much simplified (Eq. (28)), becoming 
independent of the transmission coefficients except at the basic 
frequency. A simple rule can then be used to obtain the mean- 
square pressure in retarded coordinates for the transient 
source sequence from that in source coordinates for a single 
permanent source: Make a formal transformation of the latter 

to retarded coordinates and multiply by jjL + MX2/K.J. On exam- 
ining Lighthill's work we found that his original result had 
been corrected by Ff owcs-Williams by applying a multiplicative 
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factor of [l + MX 2 /R}. If . we assume that this corresponds to 
our transient source effect, then the Lighthill - Ffowcs-Williams 
result would be applicable in the limit when each transient 
quadrupole goes through many cycles. 

It can readily be shown that, for a pulsating source at 
the center of the jet, the distance traveled downstream from 
the nozzle during one cycle is equal to the jet thickness divid- 
ed by the Strouhal number. Thus if the Strouhal number, S^, is 
0.1 the source travels ten jet diameters in completing one 
cycle. For S^, = 0.2 the source completes a cycle while moving 
downstream five jet diameters. Thus for low S-p (or low fre- 
quency) it is difficult to justify the assumption that the 
transient source passes through many cycles. However, we have 
an equation (Eq. (27)) which provides for transient sources 
passing through a reduced number of cycles down to one cycle. 

In this equation is a factor which becomes a delta-function 
in the limit as n — *■ 00 : 


2 


(w/u> 0 ) 2 sin 2 (2irnw/« 0 ) 
tr^l n - 1)2 


(109) 


In Fig. 23 we examine this factor as a function of w/cjq for 
two cycles (n = 1) and for one cycle (n = -§-). (2n must be an 

integer.) Also, we examine a related factor,, not derived here, 
which applies to half -cycles: 


8 (w/«q) 2 cos^ Jj^(m ww/cu q) J 

2 ( 2 / 2 2 
ir m (. w /w q - 1 ) v • 


F ’ 
8 


( 110 ) 
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I 



"/"o 


Pig. 23 PLOTS OP EQS. (109) AND (110) 
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(In Eq. (110) m must be odd, and m = 1 corresponds to one-half 
cycle.) 

Fig. 23 suggests that, if the pulsating source passes 
through only one cycle or one-half cycle in its lifetime, the 
factor deviates materially from the delta-function. In some 
such cases the delta-function may still be an adequate approx- 
imation. In others it may be necessary to use the more general 
equation (Eq. (27)) instead of the simpler one (Eq. (28)) which 
applies when n approaches infinity. There may be an alterna- 
tive way to define amplitude vs. frequency for the noise- 
generating sources, permitting use of the simpler equation in 
all cases, but that is beyond the scope of this study. 



CONCLUSIONS 


(1) The acoustical source traveling with the local fluid 
produces two types of disturbances. One we will term 
true sound, which consists of true waves carrying energy. 
The other we will call pseudo-sound (characterized by 
exponential decay in our analysis) which transmits energy 
only by interference with reflected waves. At a Mach 
number of zero all the far-field noise originates as true 
sound at the source. However at a Mach number as low as 
0.3 most of the far-field noise originates as pseudo- 
sound at the source. 

(2) A review of the methods for theoretically constructing 
quadrupoles in a shear layer suggests that there is no 
elementary reason for discarding dipole and source radia- 
tion in our analysis. If mean velocity varies across the 
jet, then the absence of any net source strength in the 
entire jet does not necessarily make source radiation 
unimportant. Similarly the absence of any net dipole 
strength does not result in purely quadrupole radiation. 

The dominance of quadrupole radiation in jet noise appears 
to be a feature of Lighthill’s simplified model of the 
jet, but it is not necessarily characteristic of real jets. 

(3) Moving a source off the centerline of the jet generally 
has little effect on the noise radiated as long as the 
source remains within the central high-velocity region of 
the jet. However, a source of the same strength halfway 
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out in the shear layer has greatly reduced noise radia- 
tion (which is probably to be expected because of the 
reduced speed of convection). 

(4) The noise radiated by sources on the centerline of cir- 
cular cylindrical and two-dimensional jets is compared for 
similar velocity profiles, Mach numbers and jet thicknesses. 
The circular jet cases show pressure lobes which are broader 
and rotated in closer to the downstream jet axis. Also, 

the variation of far-field noise/(f requency) with Strouhal 
number is greater for the circular jet. Here, source 
strength, defined as maximum volume introduced per unit 
time, is fixed. 

(5) The effect of raising the jet temperature above ambient 
with fixed jet velocity and source strength is illustrated. 
The Mach number based on ambient temperature is 0.7 for 
the examples given, and Strouhal numbers up to 0.2 are 
considered. The major effects are a reduction in magni- 
tude of the mean-square pressure, and rotation of the 
pressure lobes away from the downstream jet axis. 

(6) A primary effect of changing from a source to a quadrupole 
(both on the centerline) is that pressure peaks generally 
rotate in towards the downstream jet axis (see Appendix). 
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Appendix: DIPOLES AND QUADRUPQLES IN A TWO-DIMENSIONAL JET 

Equations for Mean-Square Pressure 

In our previous reports2”4 the acoustical source (i.e. 
monopole source) on the jet centerline has been taken as the 
noise generator. Since quadrupoles are frequently considered 
as the primary noise generators in a jet, we have extended the 
analysis to include dipoles (three types) and quadrupoles (six 
types), all located on the centerline of a two-dimensional jet. 
Only the results of the analysis will be given here. Details 
of the derivation are omitted since they follow the methods 
described earlier (see Section 2), methods which were devel- 
oped for the source in Refs. 2-4 and which are applied in 
detail for the source in the shear layer in Section 4 of this 
report. Following are the expressions for far-field mean-square 
pressure (Ap^) for these singularities, including the source. 

The source strength is here taken as the maximum volume 
per unit time introduced by the source. (Lighthill uses rate of 
change of mass flow as source strength.) A source strength fac- 
tor A is defined such that source strength = 8ff^A. The corres- 
ponding dipole strength factor is then A d = A^ where i is the 

2 

distance separating a source and a sink, each of strength 8v A. 

2 “T 

' Similarly, the quadrupole strength factor is A q = \l . Ap s 2 is 
the mean-square pressure produced by a source of unit strength 
factor. 
source 

o 2 2 

ip 2 = A 2 Ap 2 = A 2 - — — (z/R) 2 [s" 2 («)<-S' 2 («)] (A-l) 

S RZjl.MxVnl 5 
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x-axis dipole 

Ap 2 = Ap s 2 A d 2 w ’ 2 (X2/R) 2 /(1+MX 2 /R) 2 
y-axis dipole 

Ap 2 = Ap s 2 A d 2 «* 2 (y/R) 2 /(l+Mx 2 /R) 2 


(A-2) 


(A-3) 


z-axis dipole 


ap 2 = Ap s « d 


5.2.2 V' 2 <“>*V 2 <“> 

w f 


S"*(«0+S ,2 (iO 


1 - 


(x 7 /R) 2 +(y/R) 2 


(1+Mx 2 /R)‘ 


x-x quadropole 

Ap 2 =Ap s 2 A q 2 «' 4 (x^R^/U+Mx^R)* 
y-y quadrupole 

Ap 2 = Ap s 2 A 2 w ' 4 (y/ R) 4 / ( I+MX 2 / R) 4 


(A-4) 


(A-5) 


(A-6) 


z-z quadrupole 


Ap 2 = Ap s 2 A q 2 w ' 4 


1 - <x 2 /R) 2 -»-(y/R) 2 ^ 


(1+Mx 2 /R) 2 J 


(A-7) 


x-y quadrupole 


Ap 2 = Ap g 2 A q 2 w* 4 (x 2 /R) 2 (y/R) 2 /(l+Mx 2 /R) 4 


(A-8) 


x-z quadrupole 

^ . ^5 A 2 r .« y 2 <->*V 8 <«> 

S Q S" 2 (w)+S ,2 (u,) (1+MX2/R) 2 


1 — 


(x 2 /R) 2 +(y/R) 2 
(1+MxVR) 2 


(A-9) 
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y-z quadrupole 

(x /R) 2 +(y/R) 2 
1 

(1+Mx^/R) 2 
(A-10) 

Here, p is the mass density of the ambient air, w is the gener- 
ating frequency of the transient singularities, and <*>' is u/c 
where c is the velocity of sound. Xg, y and z are respectively 
the streamwise, lateral and normal distances of the observation 
point from the noise-generating region of the plane jet, and 
R = Vx 0 2+y 2 +z2\ yj i s the Mach number of the transient sources 
(or dipoles or quadrupoles) relative to the ambient air. . 

S' and S" are related to the in-phase and out-of-phase 
transmission coefficients for symmetrical disturbances of the 

jet (zero normal displacement at the jet centerline). S ' and 

■ A 

S^" are similarly related to transmission coefficients for anti- 

symmetrical disturbances (zero pressure at the jet centerline). 

For all of these S-terms programmed computations are desirable. 

(Expressions for S'.S" were derived in Ref. 2. These and S. ’ , 

. A 

S A " follow as special Cases (i.e. for z g = 0) of coefficients 
derived in Section 4. Using the notation of Section 4, Eqs. 

(29), (30) and (32) in particular, we find that 


•>2 _ 


Ap = Ap 


o ^ , S ,, 2 (w) + S ,2 («) 

2 A 2 ..,4 A A 


(y/R)‘ 


S A q w 


S" 2 ( w)+S* 2 ( w ) ( 1+Mx /R) 2 

2 


, _ + ^l-MsKj) 2 ^ 2 [sBO^sM i/Ul>L ] g -g(Ci)L jfj 


Den. 


S " = 

V 


SBd-M^) 




(A-ll) 


Den. 
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where the u 


*f ’ 


values of Eq. (22) are needed here. 


These also enter the denominator term. Den., Eq. (25).) 


Discussion of Equations 

Of the nine additional noise generators considered here, 
five (those involving only x and y axes) can be obtained by 
differentiation of the source far field. Three of the remain- 
ing four require anti-symmetric transmission coefficients, and 
the fourth (the z-z quadrupole) retains symmetric transmission 
coefficients but is not obtained by simple differentiation in 
the far field. 

Since for the two-dimensional jet we are examining pri- 
marily the Ap^ in the plane y = 0, we are not now concerned with 
singularities having a y-axis. Also, we are more interested in 
quadrupoles than in dipoles, and this leaves the x-x, z-z and 
x-z quadrupoles for our attention. For simplicity, we consider 
now the x-x and z-z longitudinal quadrupoles and show in Fig. 

24 the sum of the Ap 2 values for these two (Eq. (A-5) plus Eq. 
(A-7) ) . (If they wet re assumed to pulsate in phase a source-like 
Ap^ would result.) 

Fig. 25 shows for comparison purposes the Ap 2 for a 
source (Eq. (A-l)). 


Discussion of Figures 

These figures show radial plots of Ap 2 for jet Mach num- 
bers 0.3 and 0.7, for a source and for the j_(x-x) + (z-z)] quad- 
rupole (both located on the jet centerline in a jet with linear 
shear layers extending in to the centerline). The Strouhal 
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i 

A 2 n 2 / ( 2. 2 ,4 2 2. 
A p R /(2ir A q <u' p w ) 


Fig. 24 MEAN-SQUARE PRESSURE IN THE FAR FIELD (POLAR PLOT): 

' 5 — 5 

Ap (x-x quadrupole) + Ap (z-z quadrupole). 
Quadrupole at centerline; maximum jet Mach numbers 


0.3 and 0.7 
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number, S T , is defined as 

S T = (w/2w)(jet thickness)/( jet velocity) 

where u» is the source or quadrupole frequency. Results are 
shown for S^, = 0.1, 0.2, 0.5 and 2.0. The radial line with 
e on one side and w on the other indicates the boundary be- 
tween the region in which the disturbances leaving the singu- 
larity might be termed "pseudo-sound" and are exponentially 
decaying in z (e) and the region in which disturbances leaving 
the singularity are "true sound", or true waves (w). At a 
Mach number of zero all far-field noise is produced by dis- 
turbances which leave the noise generator as true sound. How- 
ever, as noted in other examples, at a Mach number of only 0.3 
most of the far-field noise is produced by disturbances ori- 
ginating as pseudo-sound. 

—X* 

The amplitudes for Ap^ are of interest (at present) only 
for comparing like Strouhal numbers at different Mach numbers 
for a given type of noise generator. It can be seen that for 
a fixed strength of the noise generator Ap2 rises rapidly with 
increasing Mach number. Before attempting to construct a power 
law for noise level as a function of Mach number other factors 
must be considered. For example, the dependence of noise gen- 
erator strength on Mach number, the possibility that noise 
generator type depends on Mach number, the conversion of gener- 
ating frequency to observed frequency, the contributions of 
various portions of the jet (which travel at different velo- 
cities) are some of these factors. - 

A primary effect of changing from sources to quadrupoles 
is that the pressure peaks generally rotate in towards the 


downstream axis 
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